Measurement of the ZZ production cross section in pp collisions at yfs =1.96 TeV 
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We present a new measurement of the production cross section o(pp — > ZZ) at a center-of-mass 
energy yfs = 1.96 TeV, obtained from the analysis of the four charged lepton final state t^l~ I ~ 
{£, i — e or fi). We observe ten candidate events with an expected background of 0.37±0.13 events. 
The measured cross section cr(pp — >■ ZZ) = 1.261'Iq'37 (stat) ± 0.14 (syst) pb is in agreement with 
NLO QCD predictions. This result is combined with a previous result from the ZZ — ^ t^l~vv 
channel resulting in a combined cross section of cr(pp ZZ) = lAOto^^ (stat) ± 0.14 (syst) pb. 

PACS numbers: 12.15.Ji, 13.85.Qk, 14.70.Hp 



Studies of the pair production of electroweak gauge 
bosons provide an important test of electroweak theory 
predictions. The production of pairs of Z/^* bosons 
has the smallest cross sections for any standard model 
(SM) diboson process not involving the Higgs boson. 
The next-to- leading order (NLO) SM prediction for the 
Z/^*Z/^* production cross section in pp collisions at 
the Fermilab Tevatron Collider at ^/s = 1.96 TeV is 
a{pp Z/-f*Z/j*) = 1.4± 0.1 pb [1]. This cross section 
is evaluated in a high mass region where the masses of 
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(Z/7*)i and (Z/7*)2 are greater than 70 GeV and 50 
GeV, respectively. A correction factor of 0.93, derived 
using PYTHIA [ij, is used to convert the measured cross 
section for Z/j*Z/'j* into that for ZZ production. Stud- 
ies of this process are important not only to further test 
the SM, but also for Higgs boson searches. Specifically, 
if the Higgs boson has a mass greater than 180 GeV, it 
will have a significant branching fraction into Z boson 
pairs. Thus, in that context, SM Z/j*Z/j* production 
is an important background to Higgs boson searches. Be- 
yond the Higgs sector, the observation of an unexpectedly 
high cross section could indicate the presence of anoma- 
lous ZZZ or ZZj couplings 3J or the existence of extra 
dimensions ^ or exotic particles. 

Previous investigations of Z/j*Z/j* production have 
been performed both at the Fermilab Tevatron pp and 
the CERN e+e" (LEP) Colliders 0. The CDF col- 
laboration reported evidence of ZZ production with a 
significance of 4.4 standard deviations from combined 
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ZZ -)> and ZZ l+l'vv searches and mea- 

sured a production cross section of a(ZZ) = 1.4^0.6 pb 
with 1.9 fb~^ of integrated luminosity [6|. The DO collab- 
oration reported an observation of ZZ — )■ ~ {t, 
I = e or /i) with 1.7 fb^^ of data and measured the pro- 
duction cross section to be (t{ZZ) = 1.75j;j!);gg (stat) ± 
0.13 (syst) pb [7|. That result was combined with a 
previous ZZ — >■ M analysis [Si] and an analysis in the 
ZZ t+t-vv channel giving a cross section of 
a{ZZ) = 1.60 ± 0.63 (stat)^Q (syst) pb with a sig- 
nificance of 5.7 standard deviations [TJ. 

In this Article, we present a measurement of Z/^* 
boson pair production with subsequent decays to either 
electron or muon pairs, resulting in final states consisting 
of four electrons (4e), four muons (4/z), or two muons and 
two electrons (2^2e) [13] • We accept events which have 
more than four leptons, however we only use the four 
leptons with highest transverse momenta in construct- 
ing kinematic variables. Compared to previous publica- 
tions d, 01 we use a larger dataset and more inclusive 
selection criteria to achieve a reduction of a factor of 2.5 
for the statistical uncertainty which dominates the exper- 
imental cross section determination. The higher statistics 
opens the possibility of studies of Z/^*Z/^* production 
properties, and we present for the first time differential 
distributions for the final state particles. Data used in 
this analysis were collected with the DO detector at the 
Fermilab Tevatron pp Collider at -/s = 1.96 TeV be- 
tween April 2002 and March 2010 and correspond to an 
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integrated luminosity of 6.4 ± 0.4 fb" 

The DO detector [l2| consists of a central tracking sys- 
tem, a calorimeter, and a muon detection system. A 
silicon microstrip tracker (SMT) and a scintillating fiber 
tracker (CFT) comprise the tracking system, which pro- 
vides coverage for pseudorapidity |77det| < 3 [l^]- The 
tracking systems are located within a 2 T superconduct- 
ing solenoidal magnet. Located immediately before the 
inner layer of the calorimeter is the central preshower de- 
tector (CPS), consisting of approximately one radiation 
length of absorber followed by three layers of scintillat- 
ing strips. Calorimetry is provided by three liquid argon 
and uranium calorimeters. The central calorimeter (CC) 
provides coverage for |?7dct| < 1-1, while the two end- 
cap calorimeters (EC) extend coverage to |?7det| < 3.2. 
The calorimeters are sectioned in order of increasing dis- 
tance from the collision point. The section closest to 
the collision region is the electromagnetic section (EM), 
while farther away are the fine hadronic (FH), and the 
coarse hadronic (CH) sections. A muon system surrounds 
the calorimeters, consisting of three layers of scintillators 
and drift tubes and 1.8 T iron toroidal magnets, covering 
hdoti < 2. 

All events used in this analysis are recorded after sat- 
isfying a mixture of single and dilepton triggers. Due to 
the high transverse momentum of the Z /^* decay prod- 
ucts and the number of leptons in the final state, the 



trigger efficiency exceeds 99%. 

The 4e channel requires the presence of four electrons 
with transverse energies Et > 30, 25, 15, and 15 GeV, 
respectively. Electrons can be reconstructed in either 
the CC region or in the EC region, however at least two 
electrons must be in the CC region. Electrons must be 
isolated from other energy clusters in the calorimeter and 
have a large fraction of their energy deposited in the EM 
section of the calorimeter. Electrons in the CC are re- 
quired to satisfy identification criteria based on multivari- 
ate discriminants which use calorimeter shower shape, 
CPS, and tracking information. Several of these parame- 
ters are inputs to a neural network (NN) , which is used to 
enhance electron purity. Electrons in the CC are required 
to have a matched track in the central tracking system. 
Electrons in the EC are not required to have a track 
matched to them due to deteriorating tracking coverage 
for |77dot| > 2, but must satisfy additional shower shape 
requirements as well as pass tighter NN selections. With 
no requirement applied on the charge of the electrons 
to increase selection efficiency, three possible Z/'y*Z/^* 
combinations can be formed for each 4e event. Only 
events having an invariant mass pair > 70 GeV and the 
other pair > 50 GeV are considered. Finally, events 
are split into three categories, depending on the num- 
ber of electrons in the CC region. Subsamples with two, 
three, and four electrons in the CC are denoted as 4e2C, 
4e3c, and 4e4c, respectively. This splitting is performed 
because these subsamples have different levels of back- 
ground contamination. 

For the 4^ channel, muons are identified as track seg- 
ments in the muon detector matched to a central track or 
as a central track matched to a pattern of calorimeter ac- 
tivity consistent with passage of a high momentum muon. 
Muons identified in the muon system must satisfy quality 
criteria based on scintillator and wire information, and 
be synchronous with the beam crossing time to reject 
background from cosmic rays. At least three muons in 
the event must be isolated. Muon isolation is dependent 
upon two cone-based variables. The first variable, Tnaio, 
is the sum of the transverse momentum associated with 
tracks in a cone of radius ATZ — \J {A-qY + (A^)^ — 0.4 
centered on the muon track. The second variable, Cnaio, 
is the transverse energy measured in the calorimeter, in 
an annulus between ATZ = 0.1 and ATZ = 0.4 centered 
on the muon track. Muons with muon system recon- 
structed tracks are considered isolated if Tnaio is less 
than 4 GeV. For muons with only a calorimeter signal 
or where the muon system provides track segments only, 
a tighter isolation requirement is used: THaio/Pj- < 0.09 
and (Cnaio — 0.005£)/p^ < 0.09, where is the trans- 
verse momentum of the muon track, and C represents 
the instantaneous luminosity (in units of 10"^° cm~^s~"'^, 
£ can reach «300) which is introduced to account for 
the occupancy increase due to multiple pp interactions 
at higher luminosities. We require that the four most en- 
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ergetic muons have ordered transverse momenta > 30, 
25, 15, and 15 GeV, respectively. The difference between 
the distances of closest approach (dca) to the pp inter- 
action point along the beam axis for any pair of muon 
tracks are required to be < 3.0 cm. The three possible 
Z/^*Z/"f* combinations per event formed without con- 
sidering muon charge are considered. Candidate events 
are selected when at least one of the three possible com- 
binations satisfies the same dilepton invariant mass re- 
quirements applied in the 4e channel. 

For the 2/i2e channel, one electron and one muon must 
have Et{pt) > 20 GeV, while the other two leptons must 
have Et{pt) > 15 GeV. All muons and electrons must 
satisfy the lepton selection criteria defined for the 4e and 
4/z final states, except that only one muon must satisfy 
the isolation requirements imposed in the 4/z final state. 
In addition, electrons and muons are required to be spa- 
tially separated by ATZ > 0.2. This requirement is ap- 
plied to remove Z — >■ ^/i background where the muons 
radiate photons leading to events with two muons and 
two trackless electron candidates. Events from this chan- 
nel assume that the muon pair originated from one Z/j* 
and the electron pair originated for the other Z/j*. The 
two same-flavor lepton pairs are required to satisfy the 
same invariant mass requirements as for the 4e channel. 
Finally, events are split into three categories depending 
on the number of electrons in the CC region. Subsam- 
ples with zero, one, and two or more electrons in the CC 
are denoted as 2/i2eoc, 2^2eic', and 2/i2e2Cj respectively. 
As in the 4e channel, this splitting is performed because 
these subsamples have different levels of background con- 
tamination. 

A Monte Carlo (MC) simulation is used to determine 
signal acceptances, efficiencies as well as the expected 
number of signal events in each subchannel. All sig- 
nal acceptances and efficiencies are evaluated after the 
high mass (> 70 GeV and > 50 GeV) requirements have 
been applied at the MC generator level. The contribu- 
tion from Z/7*Z/7* events with at least one Z/j* boson 
decaying into tau pairs is included in the signal. Events 
are generated using pythia and passed through a de- 
tailed GEANT-based 1J| simulation of the detector re- 
sponse. Differences between MC and data reconstruction 
and identification efficiencies for electrons and muons are 
corrected using efficiencies derived from large data sam- 
ples of inclusive Z ^ ££ events. 

Backgrounds to the Z/j*Z/'y* signal originate from 
events with W and/or Z bosons decaying to leptons plus 
additional jets or photons and from top quark pair (ti) 
production. The jets can be misidentified as leptons or 
contain electrons or muons from in-flight decays of pions, 
kaons, or heavy-flavored hadrons. 

To estimate the background from events with misiden- 
tifled leptons, we flrst measure the probability for a jet 
to produce an electron or muon that satisfles the identi- 
fication criteria from data. We measure this probability 



in a separate dijet data sample, selected by requiring at 
least two jets with pT > 15 GeV. We require the jet with 
largest px to pass strict jet identification criteria and we 
use the second jet to measure the probability for a jet to 
be misidentified as a lepton. The two jets are required to 
be separated in azimuth by A(j) > 3.0. To suppress con- 
tamination from T4^-|-jet events, we require the missing 
transverse energy I^t < 20 GeV 15| . The lepton identi- 



fication criteria are applied to the second jet to measure 
how often a jet mimics an electron or produced a muon. 

The probability for a jet to mimic an electron, parame- 
terized in jet Ej- and rj, is approximately 4 x 10^* for the 
case of CC electrons with a matched track and approxi- 
mately 2 X 10"'^ in the case of EC electrons for which no 
track match criterion is applied. The probabilities for jets 
to be misidentified as electrons are then applied to jets in 
eee-|-jets and ^^e+jets data to determine the background 
to the 4e and 2/Lt2e channels, respectively. This method 
takes into account contributions from Z-f jets, .Z-l-7-l-jets, 
WZ+jets, WW+iets, M^-l-jets, and events with > 4 jets. 
However, it counts the contribution from Z-|-jets twice. 
A correction is measured using data, amounting to ap- 
proximately 10%. 

The probability for a 15 GeV (100 GeV) jet to produce 
a muon of > 15 GeV is approximately 7x 10~^ (10~^) 
without requiring muon isolation, and approximately 
4 X 10^^ (2 X 10^"^) when the muon is required to be 
isolated. The probabilities for jets to contain a muon are 
applied to jets in ^^+jets and ee-|-jets data to estimate 
the background for the 4/i and 2/Lt2e channels. 

The background from tt production is estimated from 
simulation with ALPGEN [16] generated events interfaced 
to PYTHIA [2] . 

Another possible background in the 4/i and 2/i2e chan- 
nels is from cosmic ray muons. The probability for cos- 
mic ray muons to cross at the interaction region near 
the time of the pp collision is small, nonetheless we esti- 
mate this background using data. The estimation is done 
by reversing combinations of the 4/i sample selection re- 
quirements, such as scintillator timing and dca criteria. 
This procedure yields rejection factors which are then 
applied to a cosmic ray enhanced data sample. The re- 
sulting background from cosmic rays in the Afj, and 2/i2e 
samples is less than 0.01 event for each channel. 

We also estimate the contribution of Z/'j*Z/j* with 
low invariant mass lepton pairs (< 70 GeV and < 
50 GeV) that pass the kinematic selection criteria due 
to detector and reconstruction effects. This migration 
contribution is found from our signal MC where we se- 
lect events that fail the generator level mass selection. 
This small contribution is corrected for in the cross sec- 
tion measurement. 

Table H] summarizes the expected signal and back- 
ground contributions to each subchannel, as well as the 
numbers of candidate events in data. The systematic un- 
certainty for the signal yield is dominated by a 6% uncer- 
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TABLE I: The expected number of Zj^* Zj'^* and background events [tt, VK/Z/7+jets, and cosmic ray contributions], and the 
number of observed candidates in the seven Z/'y* Z/^* — > I ~ subchannels. The expected number of .2/7*^/7* events 
assumes the NLO theoretical cross section of 1.4 pb. Uncertainties reflect statistical and systematic contributions added in 
quadrature. 



Subchannel 


4e2c 


4e3c 


4e4c 


4m 


2^2eoc 


2At2eic 


2^2e2c 


Zh'Zh' 


0.31 ±0.05 


0.73 ±0.12 


0.69 ±0.11 


2.57 ±0.36 


0.24 ±0.03 


1.41 ±0.18 


2.58 ±0.33 


Zj^'Zj'y* Migration 


u.uia_o 




020+" °°® 

U.UZU_o 006 


106+°"^^ 

U.iUD„o.oi5 


U.UUZ_o 001 


002+° °°^ 

U.UUZ_o 001 


o.oost":^^ 


W/Z/7+jets 


0.065 ±0.013 


0.041 ± 0.007 0.024 ± 0.007 0.035 ± 0.015 


030+° ''" 

U.UOU_o 009 


U.UO(_o 009 


0.078tHli 


Cosmics 
it 

Observed events 





1 


2 


< 0.01 

4 


< 0.001 
001 3+° °°^° 

U.UU10_o 0009 




< 0.003 

U.Uirf»_0 0069 

1 


< 0.006 
0091+° °°''i 

U.UU»1_0 0039 

2 



tainty on the luminosity measurement ll|, the theoreti- 
cal cross section uncertainty of 7%, and the uncertainty 
on the four-lepton reconstruction efficiencies of ~ 10%. 
Additional smaller systematic uncertainties arise from 
modeling energy and momentum resolutions and from 
MC modeling of the signal kinematics. A systematic un- 
certainty of 20% on the jet-to-electron misidentification 
probability is estimated by varying the selection crite- 
ria of the control samples. Systematic uncertainties on 
background from jets containing a muon arise from the 
40% uncertainty in measured misidentification rates and 
from the limited statistics of the data remaining in the 
samples after selection. The it background systematic 
uncertainty includes the 7% uncertainty on a{tt)^ as well 
as contributions from the variation in cross section and 
acceptance originating from the uncertainty on the mass 
of the top quark. 
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FIG. 1: Distribution of the dilepton masses compared to the 
expected signal and background. 

The expected number of signal and background events 
are 8.73 ± 1.22 and 0.37 ± 0.13, respectively. We observe 
a total of ten candidate events, three in the 4e subchan- 
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FIG. 2: Distribution of Z/'y*Z/'y* pr compared to the ex- 
pected signal and background. 



nel, four in the 4/i subchannel, and three in the 2/j,2e 
subchannel. 

Figures [THl] show four kinematic distributions of the 
data compared to the expected signal and background. 
In the eeee and ^/i^/i subchannels there can be up to 
three possible pairings of the four leptons which satisfy 
the invariant mass requirements used to select candidate 
events. If two or more combinations satisfy the invariant 
mass requirements we select the one in which both dilep- 
ton pairs have an invariant mass closest to the nominal Z 
boson mass for the distributions shown in Figs. [Hand |3l 
Figure [T] shows the distribution of dilepton masses (two 
entries per event) , Fig. [2] the transverse momentum of 
the Z/"f*Z/^* system. Figure [3] displays the azimuthal 
angle (/)dccay, i-e. the angle through which the lepton side 
of one of the Z/^* boson decay planes is rotated into 
the lepton side of the other Z/^* boson decay plane, as 
measured in the Z/^*Z/^* center-of-mass frame. This 
angle is discriminating against background for high mass 
Higgs bosons. The construction of ^idecay used in this 
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FIG. 3: Distribution of the azimuthal angle <^dccay for the 
decay planes of the Z/j* bosons compared to the expected 
signal and background. 
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FIG. 4: Distribution of four-lepton invariant mass compared 
to the expected signal and background. 

Article follows the definition in Figure |4] displays 

the invariant mass of the Z/j*Z/'^* system. Additional 
differential distributions and event information for the 
selected sample of events are shown in (isj . 

The distributions shown are consistent with the expec- 
tation of a Z/'-f*Z/j* signal and small background. We 
therefore proceed to measure the pp ^ Z/'j*Z/j* pro- 
duction cross section a. Using the following likelihood 
function: 



L(ivf^M,) = ^ 



N° 



13 ■ 



(1) 



where A^°'^'' is the observed number of events given an 
expected signal and background yield of 



Here, Aj is the acceptance times efficiency, Cj is the in- 
tegrated luminosity, Bj is the branching fraction, and 



— a X X Bi X d 



bkgd 



(2) 



jybkgd ^-j^g expected background for subchannel j. The 
cross section a is obtained by minimizing — ln(L). The 
statistical uncertainty on a is obtained by varying the 
— ln(L) by half a unit above the minimum. Systematic 
uncertainties are propagated to cross section uncertain- 
ties via variations in the likelihood function due to each 
independent systematic source. These likelihood varia- 
tions are then summed in quadrature to obtain the total 
systematic uncertainty. 

The production cross section is measured to be 
aipp ->■ Z/j*Z/-f*) = 1.33t^;^S5 (stat) ± 0.12 (syst) ± 
0.09 (lumi) pb. This result is consistent with the SM 
prediction of 1.4 ± 0.1 pb. The total uncertainty reflects 
an improvement by a factor of approximately 2.5 rela- 
tive to our previous four charged lepton measurement j3| . 
Based on this result we also quote a measurement of the 
on-shell (T{pp — > ZZ) cross section. Using the conversion 
factor of 0.93 found from simulation, we measure cr(pp — > 
ZZ) = l.2Atoil (stat) ± 0.11 (syst) ± 0.08 (lumi) pb. 

The significance of the observed event distribution is 
found by using a negative log-likelihood ratio (NLLR) 
test statistic defined as —2ln{Ls+B/LB), where Lb and 
Ls+B are Poisson likelihood functions for background 
and signal plus background, respectively . As in- 
put we use the expected numbers of events from sig- 
nal and background, separated into the seven subchan- 
nels, compared to the observed numbers of data events. 
The significance is obtained by generating many pseudo- 
experiments which are created by varying the signal and 
background around their central predicted values, thus 
creating a distribution of NLLRs. The mean numbers 
of expected signal and background events per pseudo- 
experiment are varied according to their systematic un- 
certainties. The method gives the probability (p-value) 
of the background fluctuating to give the observed yields 
or higher. In 2x10^ background pseudo-experiments, we 
find zero trials with an NLLR value smaller or equal to 
that observed in data. This gives a p-value of less than 
10~^. The equivalent probability for a Gaussian distri- 
bution is greater than 6 standard deviations. 

Finally, this result is combined with the result from 
the independent ZZ — )■ i^i~vv analysis The com- 
bination is done by adding the ZZ — i^£~vv results in 
dielectron and dimuon final states to our likelihood calcu- 
lation as additional channels. Correlations of systematic 
uncertainties are accounted for between the two analyses. 
The combined result is a{pp ZZ) = lAOta^tr (stat) ± 
0.14 (syst) pb. 

In summary, the Z/j*Z/"f* cross section in pp interac- 
tions at ^5=1. 96 TeV is measured to be 1.331q 4q (stat)± 
0.12 (syst) ± 0.09 (lumi) pb. The on-sheU ZZ pro- 
duction cross section is 1.24lo;37 (stat) ± 0.11 (syst) ± 
0.08 (lumi) pb. The new DO combined result is (j{pp — )■ 
ZZ) = 1.401^-3? (stat) ± 0.14 (syst) pb. These re- 
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suits constitute the most precise measurement to date 
of the pp — )■ Z/^*Z/^* and pp — > ZZ cross sections 
and demonstrate sufficient statistics for an examination 
of Z/^*Z/^f* kinematic distributions. The kinematic dis- 
tributions of the 10 observed events are consistent with 
the SM predictions. 
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SUPPLEMENTAL MATERIAL 

Figures [5]-[9] show various kinematic quantities for data and for the expected signal and background. Figures [5] and[6] 
show individual lepton pT and r^dot, where the leptons are ranked from highest transverse momentum to lowest [13|. 
Figure [7] shows A(j> and ATZ between two leptons (two entries per event). Figure [8] shows the dilepton pair pr ranked 
from highest transverse momentum to lowest. In the eeee and ^/i^/i subchannels there can be up to three possible 
pairings of the four leptons which satisfy the invariant mass requirements used to select candidate events. If two or 
more combinations satisfy the invariant mass requirements we select the one in which both dilepton pairs have an 
invariant mass closest to the nominal Z boson mass. For the 2/z2e subchannel, we use the only valid combination, 
where one Z/'j* decays to the ee pair and the other Z/j* decays to the fifi pair. Figure [9] shows the distribution 
of dilepton mass M2 versus dilepton mass Mi, where Mi is the dilepton mass associated with the lepton pair that 
has the higher px- Figure [TU] gives an illustrative description of the angle 4>decay Tables HfflXIl give various measured 
quantities of the ten candidate events. 




50 100 150 200 50 100 150 200 
Lepton 3 p.^ (GeV) Lepton 4 p.^ (GeV) 



FIG. 5: Distributions of transverse momentum compared to the expected signal and background for the (a) leading, (b) second, 
(c) third, and (d) fourth leptons. 
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FIG. 6: Distributions of ?7det compared to the expected signal and background for the (a) leading, (b) second, (c) third, and 
(d) fourth leptons. 




FIG. 7: Distributions of and AT?, between leptons compaxed to the expected signal and background. For the eeee and /ti/U/ti/U 
channels, the combinations shown are those which are most consistent with a Z boson mass hypothesis of 91.2 GeV. 
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FIG. 8: Distributions of Zj'^* pr for lepton pairings with (a) highest and (b) lowest Z/'^* pr compared to the expected signal 
and background. 
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FIG. 9: Dilepton mass pairing compared to the expected signal and background. The lines indicate where the invariant mass 
requirements are applied. 




FIG. 10: The definition of the azimuthal angle </>dccay = 4> £ [0,27r]. l\ and I2 represent the anti-leptons from Z bosons 1 and 
2, respectively. Obtained from [ITj]. 
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TABLE II: Summary of the properties of the 1st eeee candidate event. 



Run Number 231347 
Evout Number 25076242 





('i 


<>> 




(-4 


PT (GeV) 


107 ±4 


59 ±3 


52 ±2 


17 ± 1 




0.66 


0.25 


-0.64 


-0.85 


4> (radians) 


4.09 


1.08 


0.46 


2.62 


Zvtx (cm) 


3.4 


3.4 


3.4 


3.4 


Charge 


+1 


+1 


-1 


-1 










Dilepton mass (GeV) 


89 


± :^ 


60 


± 2 


Z/V'pj (Go\-)' 


iiO 


^ 1 


iOG 


^ 3 


4-lepton mass (GeV) 
Z/-i*ZI'y''pT (GeV) 
(GeV) 




274 ± 6 
5 ± 5 
3 ± 2 





TABLE III: Summary of the properties of the 2nd eeee candidate event. 



Run Number 
Event Number 




223736 
14448774 






ei 


£2 


63 




PT (GcV) 


83 ±3 


75 ±3 


35 ±2 


27 ± 2 


n 


0.64 


0.39 


0.85 


1.18 


(j) (radians) 


6.16 


3.80 


3.83 


1.40 


«vtx (cm) 


-19.2 


-19.2 


-19.2 


-19.2 


Charge 


+1 


+1 


-1 


-1 




et 


H 


4 




Dilepton mass (GcV) 


98 


± 3 


88 


± 4 


Zh-pT (GeV) 


63 


± 3 


57 


± 3 


4-lcpton mass (GeV) 
Zh''Zh*pT (GeV) 
(GeV) 




216 ± 5 
51 ± 3 
18 ± 13 





TABLE IV: Summary of the properties of the 3rd eeee candidate event. Electron e2 does not have a matched track, therefore 
no charge measurement. 



Run Number 
Event Number 


234030 
25104195 


ei e2 63 e4 


PT (GeV) 
f] 

4> (radians) 
2vtx (cm) 
Charge 


69 ± 3 64 ± 2 30 ± 2 27 ± 2 
-0.96 -2.00 -1.19 0.06 
5.41 2.54 0.05 2.53 
31.8 31.8 31.8 31.8 
-1 - +1 +1 




Dilepton mass (GeV) 
Z/7>T (GeV) 


98 ± 3 90 ± 3 
44 ± 3 45 ± 2 


4-lepton mass (GeV) 
ZI-i*ZI'~i*pT (GeV) 
(GeV) 


239 ± 5 

1 ± 4 

2 ± 1 
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TABLE V: Summary of the properties of the 1st mmMM candidate event. 



Run Number 246915 
Event Number 2000.3687 





/'i /'2 


/':', 


PT (GeV) 


5611" 491^ 


4112 3812 




0.09 0.06 


-0.99 0.81 


<p (radians) 


3.64 0.19 


2.05 5.96 


Zvtvi (cm) 


15.2 15.2 


15.2 15.2 


Charge 


+1 -1 


+1 -1 








Dilepton mass (GeV) 


9111° 


871? 


Z/7*pT (GeV) 


4111 


5412 


4-lepton mass (GeV) 


218.5 


+16 
-12 


Zh''Zh*pT (GeV) 


17 


+ 13 

-10 


^T (GeV) 


8 


± 6 



TABLE VI: Summary of the properties of the 2nd jijijin candidate event. 



Run Number 
Event Number 


248990 
47671351 


Ml M2 M3 M4 


PT (GeV) 
V 

<j) (radians) 
Zvtx (cm) 
Charge 


531? 461^ 451^ 451° 
-1.12 -1.96 -1.14 -1.78 
4.64 2.15 2.41 5.98 
22.4 22.4 22.4 22.4 
-1 -1 -hi -hi 


Ml Ms A»2 m1 


Dilepton mass (GeV) 


88l? 86l? 


Z/j*pt (GeV) 


4412 3ll^ 


4-lepton mass (GeV) 


202 Hf 


Z/y'-Z/y'-pT (GeV) 
(GeV) 


20 1? 
7 ±5 



TABLE VII: Summary of the properties of the 3rd HHHH candidate event. 



Run Number 
Event Number 


232216 

15136574 


Ml M2 M3 M4 


PT (GeV) 
V 

4> (radians) 
«vtx (cm) 
Charge 


1161^^ 7811^ 4212 241^ 
-0.04 -1.01 0.77 -1.93 
1.69 4.26 5.29 0.35 
17.4 17.4 17.4 17.4 

-hi -1 -1 +1 


M^M3 M2 m1 


Dilepton mass (GeV) 
Zh*pT (GeV) 


148 lig 90lg^ 
80 t% 6211? 


4-lepton mass (GeV) 
Z/y*Z/y*pT (GeV) 
(GeV) 


347 li 

IS +50 
18 _29 

3 ±2 
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TABLE VIII: Summary of the properties of the 4th ii^,^,^, candidate event. 



Run Number 222870 
Event Number .■?87r)4r)12 





/'i 


/'2 


/':! 


fM 


PT (GeV) 






221^ 


20tl 




-1.02 


-1.23 


1.94 


-1.56 


cj) (radians) 


6.20 


3.57 


3.35 


1.26 


Zvtyc (cm) 


-16.1 


-16.1 


-16.1 


-16.1 


Charge 


+1 


-1 


-1 


+1 





IJ'ttJ'2 




Dilepton mass (GeV) 


90+1° 


120 1? 


Z/7*pT (GeV) 




21 tl 


4-lepton mass (GeV) 


270 


+22 
-17 


Z/'y*Z/-y*pT (GeV) 


18 


+ 13 

-9 


(GeV) 


16 


±11 



TABLE IX: Summary of the properties of the 1st 2/i2e candidate event. 



Run Number 244006 
Event Number 24854310 





ei 


62 


yui 


yU2 


PT (GeV) 


50 ±2 


41 ±2 




45tt 


V 


0.52 


-0.76 


-0.03 


-1.30 


4> (radians) 


2.22 


0.46 


3.5 


5.41 


Zvtx (cm) 


6.5 


6.5 


6.5 


6.5 


Charge 


+1 


-1 


-1 


+1 



^t^2 Ml f^t 


Dilepton mass (GeV) 
Z/7>T (GeV) 


94 ±3 102 1" 
58 ± 2 56.1 +1 


4-lepton mass (GeV) 
Z/rZ/rPT (GeV) 
(GeV) 


236 +\l 
13 1? 
9 ±6 



TABLE X: Summary of the properties of the 2nd 2/i2e candidate event. 

Run Number 244503 
Event Number 19036212 





ei 


62 


Ml 


M2 


PT (GeV) 


69 ±3 


36 ±2 




571^° 


V 


-0.86 


-0.42 


0.37 


0.15 


(j> (radians) 


1.59 


5.27 


4.76 


0.14 


2vtx (cm) 


-11.3 


-11.3 


-11.3 


-11.3 


Charge 


+1 


-1 


-1 


+1 





6+62 


Mr 


M2 


Dilepton mass (GeV) 
Z/j'-pT (GeV) 


98 ±3 
42 ±3 


86 
78 


+11 

-8 

+10 

-7 


4-lepton mass (GeV) 

Z/'y*Z/j*pT (GeV) 
^T (GeV) 


238 

76 
72 


+15 

-12 

+10 

-7 

±51 
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TABLE XI: Summary of the properties of the 3rd 2/i2e candidate event.. Electron ei does not have a matched track, therefore 
no charge measurement. 



Run Number 


208914 


Event Number 


57115998 



PT (GeV) 


56 ±2 


16 ± 1 


72+_\l 






-2.13 


0.61 


0.80 


0.25 


(j) (radians) 


5.40 


6.00 


3.00 


1.23 


^;vtx (cm) 


23.3 


23.3 


23.3 


23.3 


Charge 




+1 


-1 


+1 





(■!(■_, 


/'i 


1 


Dilcpton mass (GeV) 
Zh'pT (GoV) 


112 ±4 

70 ±2 


79 j 
76 j 


-10 
-7 
-15 
-10 


4-lepton mass (GeV) 
Z/-i*Z/')*pT (GeV) 
(GeV) 


359 
9 
8 


+30 
-21 
+ 15 
-11 

±6 





